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Resumen 


Un objetivo de los estudios hidrológicos es la estimación del volumen 


escurrido anual (VEA) de cuencas sin aforos. Datos que son básicos en 
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el diseño de embalses. En este estudio, tal problema se resuelve por 
medio del análisis de frecuencias regional, con base en el método del 
índice de escurrimientos. Este enfoque se aplicó en la Región 
Hidrológica No. 10 (Sinaloa), México, procesando 22 registros de VEA, 
con amplitudes que variaron de 24 a 57 años. Por medio de la prueba 
de discordancia se eliminaron tres registros. La homogeneidad 
hidrológica regional se verificó con los 19 registros restantes y tres 
técnicas estadísticas: Índices de estacionalidad, linealidad de los 
momentos ordinarios y regresión lineal entre momentos de 
probabilidad ponderada. Las pruebas de pérdida de homogeneidad de 
cada registro condujeron a eliminar cuatro, por mostrar tendencia 
lineal ascendente. Los 15 registros remanentes se procesaron de forma 
adimensional, al dividir sus datos entre su volumen escurrido medio 
anual (VEMA) y se concatenaron con el método de las estaciones-años, 
formando una serie de 539 datos. Para obtener la curva de crecimiento 
regional se aplicaron tres modelos probabilísticos: Pearson tipo III, 
Log-Normal y Transformación Potencial; adoptando el de menor error 
estándar de ajuste. El escalamiento de la distribución Log-Normal 
regional se logró con una relación logarítmica entre el VEMA y el área 
de cuenca de los registros procesados. Las Conclusiones detallan el 
procedimiento establecido y destacan la importancia de los resultados 
para generar secuencias sintéticas de VEA, igualmente probables de 
ocurrir. Por ello, se recomienda la aplicación sistemática del método 
del Índice de escurrimiento en otras regiones del país, para contar con 


tal método de estimación del VEA. 
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Palabras clave: volumen escurrido anual, momentos L, prueba de 
discordancia, índices de estacionalidad, homogeneidad hidrológica 
regional, tendencia lineal, error estándar de ajuste, secuendias 


sintéticas. 


Abstract 


One of the objectives of hydrological studies is the estimation of the 
annual runoff volume (VEA, for its name in Spanish) of watersheds 
without gauging. Such data is basicin the design of reservoirs. In this 
study, the estimation is provided through regional frequency analysis, 
based on the Runoff Index method. This approach was applied in the 
Hydrological Region No. 10 (Sinaloa), Mexico, for which 22 VEA records 
were processed, with amplitudes that varied from 24 to 57 years. 
Through the Discordance Test, three records were eliminated. Regional 
hydrological homogeneity was verified with the remaining 19 records 
and three statistical techniques: seasonality indices, linearity of 
ordinary moments, and linear regression between probability-weighted 
moments. The tests of loss of homogeneity of each record led to the 
elimination of four, due to their ascending linear trend. The 15 
remaining records were processed in a dimensionless way, by dividing 
theirdata over their annual mean runoff volume (VEMA) and they were 
concatenated with the stations-years method, forming a series of 539 
data. Three probabilistic models were applied to obtain the regional 
growth curve: Pearson type III, Log-Normal, and Potential 


Transformation; choosing the one with the smallest standard error of 
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fit. The scaling of the regional Log-Normal distribution was achieved 
with a logarithmic relationship between the VEMA and the basin area 
of the processed records. In the Conclusions section, the established 
procedure is detailed and the importance of the results to generate 
synthetic sequences of VEA is highlighted (which are equally probable 
to occur). Therefore, the systematic application of the Runoff Index 
procedure as a method of estimating the VEA in other regions of the 
country is recommended. 


Keywords: Annual runoff volume, L moments, discordance test, 
seasonality indices, regional hydrological homogeneity, linear trend, 
standard error of fit, synthetic sequences. 


Recibido: 13/04/2020 


Aceptado: 11/05/2021 


Introducción 


Dimensionamiento de los embalses 
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La hidrología superficial, en relación con la estimación de los recursos 
hidráulicos, intenta contestar las siguientes dos preguntas: ¿de cuánta 
agua se dispone? y ¿cómo son sus valores extremos? Ambas 
cuestiones implican la estimación de los volúmenes disponibles de agua 
y de sus magnitudes extremas (crecientes y sequías), así como su 


ocurrencia o comportamiento a lo largo del tiempo. 


Las obras hidráulicas clásicas de aprovechamiento son los 
embalses, los cuales pueden ser divididos en grandes y pequeños; los 
primeros son de propósitos múltiples y generalmente almacenan agua 
de los años húmedos y la utilizan en los años secos; los segundos 
guardan agua de las épocas de lluvias y la liberan en el siguiente 
estiaje. El diseño hidrológico de los embalses sean grandes o pequeños 
requiere de la estimación, en su sitio de proyecto, del régimen de 
escurrimientos anuales y mensuales, para su dimensionamiento. 
Además, son necesarias las crecientes de diseño para garantizar su 
seguridad hidrológica (Carr € Underhill, 1974; Zsuffa 8: Gálai, 1987; 
Skertchly-Molina, 1989). 


Debido a la baja densidad de estaciones hidrométricas, la 
estimación del régimen de escurrimientos en el sitio de proyecto del 
embalse por lo general se realizaba antes del enfoque regional, con 
base en el concepto de similitud de cuencas. Entonces, un modelo 


hidrológico del proceso lluvia-escurrimiento calibrado en una cuenca 


5 


E 2022, Instituto Mexicano de Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


o) y) Check for updates 
OPEN ACCESS q = ] 


Tecnología y 


CienciaszAgua 


aforada se aplicaba a otra similar, pero sin datos hidrométricos 
(Cannarozzo, Noto, Viola, 8 La-Loggia, 2009). 


Hoy día, el análisis de frecuencias regional permite el uso de toda 
la información hidrométrica disponible en un área geográfica a través 
de métodos o técnicas estadísticas que ayudan a definir una región o 
zona homogénea desde un punto de vista hidrológico; entonces, 
resulta factible y confiable establecer una función de distribución de 
probabilidades (FDP) común y adimensional en relación con un valor 
Índice medio o mediano. Posteriormente, en cada región homogénea, 
se obtiene una relación empírica entre el valor índice, y diversas 
características morfológicas y climáticas de la cuenca fácilmente 
calculables, haciendo posible la generación de muestras sintéticas del 
volumen escurrido anual (VEA) en cualquier sitio de la región estudiada 
(Cannarozzo et al., 2009). Blóschl, Sivapalan, Wagener, Viglione y 
Savenije (2013) indican que la hipótesis que rige este método 
establece que el volumen escurrido medio anual puede vaniar entre 
sitios de una región homogénea, dentro de la cual la FDP es idéntica y 


está definida por la curva de crecimiento (growth curve) adimensional. 


Por otra parte, Skertchly-Molina (1989) ha expuesto una 
clasificación exhaustiva de los métodos que permiten dimensionar a los 
embalses. Tales técnicas se pueden dividir en dos tipos: (1) las que 
procesan los registros históricos de escurrimiento y (2) las que 
emplean registros sintéticos igualmente probables de ocurrir. Con 
respecto a la conveniencia de unas u otras, McMahon y Mein (1986) 


señalan que el registro histórico de escurrimientos es el menos 
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probable o factible de que se repita exactamente en el futuro. Además, 
indican que contar con alguna estimación de la sensibilidad de la 
estimación de la capacidad necesaria del embalse es muy útil para 
tomar decisiones más confiables. Lo anterior es posible haciendo uso 


de series sintéticas de escurrimientos anuales (Campos-Aranda, 2011). 


Objetivo 


Corresponde en este estudio el análisis de frecuencias regional del 
volumen escurrido anual (VEA), en millones de metros cúbicos (Mm3), 
en 22 estaciones hidrométricas de la Región Hidrológica No. 10 
(Sinaloa), México. Los resultados permiten generar secuencias 
sintéticas, de cualquier tamaño, igualmente probables de ocurrir del 
VEA, en los sitios de las estaciones hidrométricas utilizadas y en 
cualquier localidad sin datos de aforo dentro de la región hidrológica 
analizada. Se utilizó el método del índice de escurrimiento, similar al 
conocido índice de crecientes. Se trabajó con 19 registros. Después de 
aplicar la prueba de discordancia y la homogeneidad hidrológica 
regional se verificó con base en tres técnicas estadísticas. En seguida, 


se definió la FDP regional a través del error estándar de ajuste. Por 
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último, se encontró una ecuación logarítmica que relaciona el volumen 
escurrido medio anual y el área de cuenca. 


Métodos y materiales 


Método del índice de escurrimientos 


Hosking y Wallis (1997) indican que las técnicas que integran el 
llamado método del índice de crecientes (index flood method) se 
pueden aplicara cualquiertipo de datos. Tal método fue propuesto por 
Dalrymple (1960) para procesar datos de crecientes, también llamadas 
avenidas máximas o riadas (de ahí su nombre), y permite estimar un 
resumen de parámetros estadísticos procedentes de la agrupación de 


varias muestras o registros diferentes (Rahman, Haddad, 8 Eslamian, 
2014; Ouarda, 2017). 

Cuando existen registros de volúmenes escurridos anuales en NS 
sitios y en cada localidad j, el número de datos es n;¡ de volúmenes v?. 


Entonces, se tiene que W(F) con O < F < 1 es la función de cuantiles 
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de la FDP de cada sitio j. La hipótesis básica del método de índice de 
escurrimientos establece que los sitios j forman una región 
homogénea, lo cual implica que las FDP de los NS sitios son idénticas, 
excepto por el factor de escalamiento, que es precisamente el índice 
de escurrimientos (IE/). La ecuación del método es: 


VÍ(E)= IE! -q(F) conj = 1,2,...,NS (1) 


siendo q(F) la curva de crecimiento regional (regional growth curve ), 
es decir, una función de cuantiles adimensional y común a todos los 
sitios; q(F) es una FDP regional o función común de los valores V? J1E), 
Por lo general, el IE/es igual al volumen escurrido medio anual (VEMA)) 


de los datos de cada sitio, pero se puede emplear la mediana (Hosking 
s: Wallis, 1997). 


En el método del IE/ se aceptan las cuatro hipótesis siguientes: 
(1) los escurrimientos de cada sitio son independientes; (2) los 
escurrimientos de sitios diferentes son independientes; (3) las FDP en 
sitios diferentes son idénticas excepto por el factor de escalamiento, y 
(4) la forma matemática de la curva regional de crecimiento está 
perfectamente especificada (Hosking 8 Wallis, 1997). 


El análisis de frecuencias regional, incluido el 14, tiene dos 
aplicaciones fundamentales dentro de la región homogénea: la primera 
se realiza en sitios con registros cortos o no confiables, y la segunda 


en localidades sin datos hidrométricos, utilizando las ecuaciones 
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empíricas obtenidas para la relación entre el VEMA! y el área de cuenca, 
u otras características morfológicas de la cuenca o del cauce principal. 


Momentos y cocientes L 


Los momentos L (A-) son combinaciones lineales de los momentos de 
probabilidad ponderada (fr) desarrollados por Greenwood, Landwehr, 
Matalas y Wallis (1979), los cuales son parámetros estadísticos de los 
datos ordenados. Los momentos L han demostrado ser un sistema 
eficiente y robusto para el ajuste de las FDP utilizadas actualmente en 
hidrología o establecidas bajo precepto. Sus ecuaciones de cálculo son 
(Hosking € Wallis, 1997; Rao € Hamed, 2000; Stedinger, 2017): 


d= Bo (2) 
A, =2*P,—Bo (3) 
A3=6:*P,— 6-B,+Bo (4) 
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A, = 20: PB, — 30: B, +12: B,— Bo (5) 


En tal sistema se definen los cocientes (Tr) de momentos L, 
comenzando con L-Cv, que es análogo a este coeficiente, y después 


los de similitud con los coeficientes de asimetría y de curtosis, los 
cuales son: 


T, =2A/My (6) 
T¿ = 2/4, (7) 
T, =2y/A, (8) 


En una muestra de tamaño n, con sus elementos X¡; expuestos en 


orden ascendente (X, < X, < *** < X, ), los estimadores insesgados de fr, 
se obtienen con la expresión general siguiente: 


1 (i-1) (1-2) ---(i—r) 
PB, ==: 


e E a conr = 0,1,...4 (9) 
Prueba de discordancia 
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En todo análisis de frecuencias regional se deben realizar, como 
mínimo, las dos verificaciones siguientes: (1) cada muestra o serie de 
datos de un sitio debe ser revisada para buscar datos erróneos, es 
decir, valores demasiado grandes o extremadamente reducidos, así 
como repetidos, los cuales se pudieron originar en la transcripción, y 
(2) se deben comparar las muestras entre ellas y con las más cercanas 
para verificar un cambio en magnitud, por ejemplo, conforme el 
tamaño de la cuenca crece o su ubicación varía de una zona a otra de 


la región analizada (Hosking 8 Wallis, 1997; Campos-Aranda, 2010) 
(Tabla 1). 


Tabla 1. Valores críticos de la discordancia (D;) según el número de 
estaciones del grupo bajo análisis (Hosking 8 Wallis, 1997). 
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Afortunadamente, los valores erróneos, los eventos dispersos 
(outliers), la tendencia, y los saltos o cambios en la media de los datos 
son reflejados en los momentos L de la muestra (ecuaciones (2) a (5)). 
Por ello, una mezcla conveniente de los cocientes L en un solo 
estadístico (D;) que mida la discordancia entre los cocientes L del sitio 
y los promedio de grupo se ha sugerido como prueba básica para 
detectar sitios que son discordantes con el grupo como un todo 
(Hosking € Wallis, 1997). 


Suponiendo que existen NS sitios en el grupo que se analiza, se 
' e ¿7T . . 
define a u; = [el el ti] como un vector que contiene los cocientes L de 


la muestra: t2, t3 y t4 de cada sitio j, definidos por las ecuaciones (6) a 
(8). El superíndice T significa transpuesto, ya que u;, es un vector 


renglón. El vector promedio (no ponderado) del grupo será (Hosking 8 
Wallis, 1997): 


ú=50, u (10) 


La matriz A de suma de cuadrados y de productos cruzados 
estará definida como: 


A= Y, (uy -1) (uy -1) (11) 


Finalmente, la medida de la discordancia de cada sitio será: 
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Entonces, cuando D; es mayor que los valores críticos (Dc), el 
sitio será discordante con el grupo. 


Índices de estacionalidad 


Zrinji y Burn (1996), Bum (1997), Cunderlik y Burn (2002), y De- 
Michele y Rosso (2002) han sugerido que la fecha de ocurrencia media 
de las crecientes y la regularidad estacional se pueden utilizar como 
medidas de similitud de la respuesta hidrológica de una cuenca, de 
manera que tales cuencas se pueden integrar en una región con 
propósitos de un análisis regional de frecuencias. Este enfoque de 
regionalización tiene la ventaja de reservar el uso de la información 
hidrométrica para verificar la homogeneidad regional. En este estudio 
del VEA se evaluó la fecha de ocurrencia media del volumen escurrido 


máximo mensual de cada año y la regularidad estacional. 


Los estadísticos direccionales datan de principios de la década de 


1970 y constituyen una herramienta simple para definir medidas de 
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similitud procedentes de las fechas de ocurrencia de los eventos 
analizados. En México fueron introducidos por Ramírez-Orozco, 
Gutiérrez-López y Ruiz-Silva (2009). 

Existen diversas convenciones o formas de usar el círculo para 
definirlos estadísticos direccionales; en adelante se usará la que tiene 
semejanza con las horas del reloj debido a que se analizan los 12 


meses del año, entonces se tiene que enero es la una; junio, las seis; 
septiembre, las nueve, y diciembre es las 12. 


Para establecer los índices de estacionalidad, se comienza por 
transformar el mes de ocurrencia de cada volumen escurrido máximo 
mensual de cada año, a ángulo de avance en grados (6,), con 30% de 
incremento en cada mes; con 15% para enero, con 75% para marzo, 


1350 para mayo, 225% para agosto, 285% para octubre, etcétera. 


Teniendo n datos del mes en que ocurrió el volumen escurrido 


máximo mensual de cada año, las coordenadas promedio del ángulo 


medio (9) de los máximos mensuales serán: 
E= a sen(8,) (13) 


y == Eh, cos(8,) (14) 
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Por lo tanto, la dirección media 6 de la ocurrencia promedio de 


los meses del escurrimiento máximo mensual en coordenadas polares 
será: 


O =fím” E) (15) 


La aplicación de la ecuación anterior se realiza obteniendo 
primero el ángulo alfa (a) tangente de x entre y, ambas con signo 
positivo; entonces si x y y son positivas O=a, si x>0 y y<0 
9 = 180%-a, siambas son negativas 9 = 180% + a y por último, six < 0 
yy>00= 360%-q. 


Se puede esperar que cuencas con valores similares del 0 
presenten semejanzas en otras características hidrológicas 
importantes. El ángulo € puede estar relacionado con el tamaño de la 


cuenca y su localización geográfica dentro de la región analizada. 


Una medida de la variabilidad de las n fechas de ocurrencia de 
los escurrimientos máximos mensuales de cada año, en relación con el 


9, se puede estimarcalculando la resultante media, cuya expresión es: 


F=.[32+y? (16) 
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El índice de estacionalidad + es una medida adimensional de la 
dispersión de los datos, toma valores entre cero y uno. Un valor 
unitario indica que todos los escurrimientos máximos mensuales 


ocurren en el mismo mes; en cambio, un valor cercano a cero implica 
gran variabilidad de ocurrencias a lo largo de todo el año. 


Prueba de la linealidad de los momentos ordinarios 


Gupta, Mesa y Dawdy (1994) encuentran que si una región es 
homogénea con respecto a la variable hidrológica X;, entonces sus 
momentos ordinarios (MOs) muestran linealidad logarítmica con las 


áreas de cuenca (4). Por lo anterior, la expresión siguiente (Ferro 8 
Porto, 2006) permite verificar tal linealidad: 


mo, = Emlsica! 


cons = 1,2,.3,4 (17) 


siendo s el orden del momento; n, el número de elementos del registro, 
y X?(4) el dato asociado con el valor del área de cuenca. Llevando a 
una gráfica semilogarítmica, los tamaños de cuenca en las abscisas en 


la escala logarítmica y en las ordenadas, en escala aritmética, los 
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logaritmos decimales de MOs, se pueden calcular las pendientes de 
cada agrupación o nube de momentos de orden s, por medio de rectas 
de regresión lineal de tipo siguiente: 


MO, =b+m:(logA) (18) 


Prueba de la regresión lineal de los B, 


Kumar, Guttarp y Foufoula-Georgiu (1994) propusieron un test que 
verifica la homogeneidad hidrológica regional, basado en una ecuación 
similar a la anterior, pero que utiliza los momentos de probabilidad 
ponderada (f,). Esta prueba fue aplicada por De-Michele y Rosso 
(2002), y Ferro y Porto (2006). Una versión asociada con el concepto 
Kumar et al. (1994) fue empleada y verificada por Varas (2000) en 
EUA, Francia, Inglaterra y Chile, para las relaciones lineales de f, de 
un cierto orden y el orden inmediato inferior. Lo anterior significa que 
se cumple la siguiente regresión lineal en una región homogénea 


hidrológicamente: 


B, =b+m:B,_, conr = 1,2,...,4 (19) 
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en la cual b y m son la ordenada al origen y la pendiente de la recta de 
regresión lineal (ecuaciones (22) y (23)), cuya calidad del ajuste la 
define el coeficiente de correlación lineal (Ty): 


E (BB 1) (BP) (20) 


Fey = : 1/2 
ARCA AT 


xy 


Recta de regresión para una tendencia lineal 


Se considera que la variable dependiente (y) son los volúmenes 
escurridos anuales VEA; en Mm3, y los tiempos o años T; son las 
abscisas (x), en este caso iguales al ¡-ésimo valor ¡. Para probar si la 
pendiente (m) de la recta de regresión ajustada por mínimos cuadrados 
de los residuos es estadísticamente diferente de cero, se usa una 
prueba basada en la distribución t de Student, definida por las 
ecuaciones siguientes (Ostle € Mensing, 1975): 
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b=VEA=m:-T (22) 


Me a 2-qp2 (23) 
t== (24) 
Sn 
$ E (25) 
m y? (71, DM? 
n 2 
Ss? E 2 ¡1 (VEA¡-VEA¡) (26) 


(n-2) 


VEA y T son las medias; en la ecuación anterior, VEA, es el valor 
estimado con la recta de regresión (ecuación (21)); Si y Si son las 
varianzas de los errores y de la pendiente. Si el valor calculado t es 
mayor que el crítico (tc) obtenido para la distribución t de Student con 
n — 2 grados de libertad y a = 5% en una prueba de dos colas, la 
pendiente m es significativa, es decir, existe tendencia lineal. El 
problema de esta prueba es que no distingue entre persistencia y 
tendencia (Adeloye € Montaseri, 2002). Para estimar el valor de t. se 
utiliza el algoritmo propuesto por Zelen y Severo (1972), con Z = 
1.95996 para una confiabilidad (1 — a) del 95 %: 
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t. =Z+G1/v+G2/v? + G3/v? + G4/v* (27) 
donde: 


G1= (23 +72)/4 
G2= (5Z* + 162? + 3Z/96) 
G3 = (327 + 197? + 172? — 15Z)/384 


G4= (7929 + 77627 + 148227? — 19207? — 9457)/92160 


Información hidrométrica disponible 


En la Región Hidrológica No. 10 (Sinaloa), México, existen 42 
estaciones hidrométricas con registro confiable y sólo 22 tienen su 
régimen de escurrimientos no afectado por los embalses y registros 
mayores de 20 años, por ello fueron las utilizadas. En el sistema 
BANDAS (IMTA, 2002) se obtuvo toda la información de volúmenes 
escurridos anuales de tales estaciones de aforos, cuyas características 


generales se presentan en la Tabla 2, expuestas en orden decreciente 
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de tamaño de cuenca. Esta información procede del CD-1, está 
expuesta a nivel mensual y se procesará, como ya se indicó, en 
millones de metros cúbicos (Mm3). Lógicamente, los registros de las 
estaciones hidrométricas Guamúchil y Huites abarcan hasta antes del 
año en que la construcción del embalse respectivo (Eustaquio Buelna 


y Luis Donaldo Colosio) afectó su régimen. 


Tabla 2. Características generales de las 22 estaciones hidrométricas 
procesadas de la Región Hidrológica No. 10 (Sinaloa), México. 


E 
lead 
sl 
La Huerta *1970-1999 (30) 
Toahayana R. Petatlán 1958-1985 (28) 


! E 
A 
a 
1 
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R. Elota 1955-1994 (40) 1 884 417.784 


| | 
A A 
Bamícori *1951-1983 (33) 


Descripción de acrónimos: 


RP = años del registro procesado. 

NA = número de años. 

A = área de cuenca. 

VEMA = volumen escurrido medio anual. 


*registros continuos. 


En la cuarta columna de la Tabla 2 se han indicado con un 
asterisco los 13 registros que son continuos, es decir, que no presentan 
dudas para la definición de su lapso de años completos. En los nueve 
restantes, que son: Santa Cruz, Palo Dulce, Ixpalino, Toahayana, 
Tamazula, Acatitán, Choix, El Quelite y Zopilote, se tienen dos o más 
años faltantes, porlo general después del año de 1980, y por ello existe 
dificultad para definir el periodo de años por procesar. 


23 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


a Check for u dates | 
OPEN Quecess | Mm P J 


Tecnología y 


CienciaszAgua 


En la Figura 1 se muestra la localización geográfica de las 22 
cuencas que definen las estaciones hidrométricas de la Región 


Hidrológica No. 10 (Sinaloa), México, citadas en la Tabla 2, las cuales 
serán procesadas en su registro de VEA. 
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Figura 1. Localización en la Región Hidrológica No. 10 (Sinaloa), 
México, de las cuencas de las 22 estaciones hidrométricas 


procesadas. 
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Lapsos de los registros por procesar 


En la definición de los lapsos de registro de volumen escurrido anual 
por procesar se pueden establecer dos enfoques de trabajo: (1) como 
registros de datos independientes e idénticamente distribuidos, y (2) 
como serie cronológica. El primer formato obedece al requisito del 
análisis de frecuencias y permite seleccionar registros interrumpidos 
por uno o varios años. El segundo formato implica conservar las 
secuencias de años húmedos y secos, y es como se deben estudiar los 
escurrimientos; por ello, se adoptó. En la cuarta columna de la Tabla 2 
se indican los lapsos de registro continuo que serán procesados. 


Definidos tales lapsos continuos se calcularon, como se detalla 
en el inciso siguiente, los datos mensuales faltantes de los años 
incompletos, así como los años faltantes, para poder definir el volumen 
escurrido medio anual (VEMA), el cual se cita en la última columna de 
la Tabla 2. 


Deducción de datos faltantes 
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Para la estimación de los volúmenes escurridos mensuales faltantes en 
los años incompletos se siguió un método simple, el cual consistió en 
promediar el mismo mes del año anterior y del posterior, para asignar 
tal valor al mes faltante. Lo anteriorse considera aceptable debido a la 
regularidad que presenta la temporada de escurrimientos año con año. 
Con tal procedimiento se completaron los registros en las estaciones 
San Francisco, Santa Cruz, Guatenipa Il, Jaina, Palo Dulce, La Huerta, 


Toahayana, Naranjo, Choix y Pericos, en uno o dos años. 


Para los años faltantes se aplicó el transporte ponderado con el 
inverso de la distancia entre estaciones hidrométricas (Campos- 
Aranda, 2008), empleando dos o más registros. Tal fue el caso del año 
1987 en la estación El Quelite, con base en Ixpalino y Acatitán; así 
como del año 1985 de Tamazula, empleando Santa Cruz y Guatenipa 
IT. 


Cuando el transporte ponderado no se pudo aplicar, el año 
faltante se dedujo con base en la moda del registro anual, estimada 
con la distribución Gamma mixta; ese fue el caso del año 1983 en 


Urique II, utilizando 26 valores anuales. 


Método regional por aplicar 
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En la cuarta columna de la Tabla 2 se destaca la gran variedad de 
fechas de inicio y terminación de los registros de volumen escurrido 
anual (VEA). Lo anterior pone en evidencia la imposibilidad de definir 
un lapso común de amplitud adecuada (> 30 años) de registros por 
procesar de VEA en la región analizada con el método del IE, 


Ante tal realidad, se optó por aplicar el método regional de las 
estaciones-años, para procesar los registros que formarán una región 


homogénea desde un enfoque hidrológico. 


Búsqueda de la FDP regional 1: Pearson tipo III 
(PT3) 


Este modelo probabilístico tiene tres parámetros: u (ubicación), OU 
(escala) y y (forma). Si y > 0, la variación de Xes¿ <X<oo, y se tiene 
la distribución PT3 más común o difundida; si y=0 ocurre la 
distribución Normal, cuyo rango de X es —o < X <oo, y cuando y < 0 se 
presenta la PT3 reversa, con un intervalo de X de —o < X < 3 (Hosking 


8 Wallis, 1997). Esta FDP no tiene solución explícita, por lo cual el valor 
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X asociado con una probabilidad de no excedencia (F) se estima a 
través del factor de frecuencia (Kf), con la ecuación general del análisis 
de frecuencias hidrológico (Chow, 1964): 


X(F)=u+K;,:0 (28) 
la expresión de Kf es (Kite, 1977): 
Fs =Z 4 (22 Dv. +52? 62)? -(22-D1 +2 yi+ ly? (29) 


en la cual y, = y/6 y Z es la desviación normal estándar función de F; 
se estima con la aproximación siguiente (Zelen 8 Severo, 1972): 


2.515517+0.802853-w+0.010328 -w? 


z=w -_— -  — — — 
1+1.432788:w+ 0.189269:'w*+0.001308:w 


(30) 
siendo: 


w= /In(1/F3) (31) 


cuando F > 0.50 se utiliza F= 1 — Fen la Ecuación (31) y al resultado 
de la Ecuación (30) se le cambia el signo. 
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El signo de 7, define el del coeficiente de asimetría (y) en la 


Ecuación (36). Si 0<|r¿| < 1/3, se hace z= 3rn(13)? para llevar a la 
expresión siguiente (Hosking 8. Wallis, 1997): 


ES 1+0.2906:z 
Dia 7+0.1882 :22+0.0442:z3 (32) 
Si 1/3 < |rz| < 1, se hace z = 1-— |13| y se aplica: 
A 72 23 
az 0.36067 :Z-0.59567:z a VA - (33) 
1-2.78861:z2+2.56096 :z*—0.77045:z 
Finalmente: 
u= Ay (34) 
o=2,'Yr:a:T()/T(a+ 1/2) (35) 
y =2/ya (36) 


Para evaluar la función Gamma se utilizó la aproximación de 
Stining (Davis, 1972): 
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Mw) =e%:-0%72-/2r1:f, (37) 
siendo: 

a 139 571 hal) 
f4= 12w 288w? 51840w?  2488320w* 


Búsqueda de la FDP regional 2: Log-Normal (LN3) 


En este modelo probabilístico sus parámetros son E (ubicación), a 
(escala) y k (forma). En el arreglo de variables de Hosking y Wallis 
(1997), el intervalo de la variable X es -o<X<8E¿+a/k cuando k > 0, 
o <X<o sik= 0 (distribución Normal), y ¿+ a/k < X < o cuando k< 
0; además, X está relacionada con la variable aleatoria Z, que tiene 
distribución Normal estándar (Ecuación (30)), definida como: 


AP)=t4 dee) cuando k + 0 (38) 
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El parámetro de forma k es sólo función de r,, con la siguiente 
relación válida para |7,| < 0.94, correspondiendo a |k] < 3: 


= 
Q 
| 


y, 20466534 -3.6544371 -13+1.839673317-0.20360244 15 
3 1-2.01582173-12+1.2420401-12-0.21741801-1% 


(39) 
Los otros parámetros son: 
E=A -(1-e72) (40) 
E A (41) 


En la ecuación anterior, d(x) es la probabilidad de no excedencia 
hasta el valor x de la distribución Normal, igual a: 


D(x) = MN $O0)dt 


(42) 
siendo p(x) la ordenada de la curva Normal en x, igual a: 
ex? 
$) = GE (43) 
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Para evaluar la Ecuación (42) se utilizó la aproximación siguiente 
(Zelen 8: Severo, 1972): 


d(x)= 1-00 (ce, d+c,: d+ -+c¿:d”) (44) 


en la cual: 


1 
— 14+0.2316419:x 


(45) 


c, = 0.31938153, c, = -0.356563782, c¿ = 1.781477937, c, = -1.821255978 y 
C¿ = 1.330274429. 


Búsqueda de la FDP regional 3: transformación 


potencial (TP2) 


Chander, Spolia y Kumar (1978) sugieren utilizarla transformación de 
Box-Cox para normalizar los datos X;¡a z;, después aplican la Ecuación 


(28) con una distribución Normal, y por último obtienen la estimación 
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buscada X(F), con la transformación inversa. Las expresiones citadas 


son: 

P_ 
z¡= “2 cuando p+0X,>0 (46) 
z¡ = 1In(X;) cuando p =0,X,>0 (47) 


La Ecuación (46) es una transformación potencial general, en la 
cual las conversiones logarítmicas, recíproca y de raíz cuadrada son 
casos especiales, ya que q es una constante que varía de -1.0 a 1.0. 
La obtención de q no es explícita, pero un incremento o decremento 
origina el mismo efecto en el coeficiente de asimetría (Cs) de los z;; 


entonces su valor óptimo será el que lleve al Cs a cero. 


Escalante-Sandoval y Reyes-Chávez (2002) indican que la 
Ecuación (46) normaliza el Cs= 0, pero no el coeficiente de curtosis 
(Ck +3). La corrección de la curtosis se logra con una segunda 
transformación, expuesta por estos autores. De acuerdo con Chander 
et al. (1978), en el análisis de frecuencias no es necesaria la segunda 
transformación. La Ecuación (28) es ahora: 


Z(F)=zm>+2Z:zS (48) 
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en la cual zm y zs son la media y desviación estándar de los datos 
transformados, y Z se estima con la Ecuación (30). La estimación 


buscada se obtiene con la expresión: 


X(F) =[p -z(F)+ 111 (49) 


Error estándar de ajuste 


De acuerdo con Chai y Draxler (2014), es el indicador más común para 
el contraste de las distribuciones de probabilidad a datos reales; se 
estableció a mediados de la década de 1970 (Kite, 1977). La expresión 
del error estándar de ajuste (FEA) es: 


_ [E%,WEA¡-VEA¡* bs 
cia Ie pe 


VEA; son los volúmenes escurridos anuales ordenados de menor 


a mayor, cuyo número es n y VEA, los volúmenes escurridos anuales 
estimados, para la probabilidad de no excedencia (F) estimada con la 


Ecuación (51) o fórmula de Cunnane, que según Stedinger (2017) 
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conduce a estimaciones insesgadas en la mayoría de las FDP utilizadas 
en hidrología o modelo probabilístico que se contrasta. Por último, np 


es el número de parámetros de ajuste de la FDP, con tres para la PT3 
y LN3 y dos para la TP2: 


_ 1-0.40 


— n+0.20 (51) 


Resultados y su análisis 


Pruebas de homogeneidad 


Para que los resultados de un análisis de frecuencias regional sean 
confiables, los registros de volumen escurrido anual (VEA) por procesar 
deben haber sido generados por un proceso aleatorio estacionario, lo 
cual implica que no haya cambiado en el tiempo. Por lo anterior, los 


registros de VEA necesitan estarintegrados con datos independientes, 
libres de componentes determinísticas. 
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Para comprobar lo anterior se aplicaron siete pruebas 
estadísticas, una general, el Test de Von Neumann y seis específicas: 
dos de persistencia (Anderson y Sneyers); dos de tendencia (Kendall y 
Spearman); una para la variabilidad (Bartlett), y la última relacionada 
con el cambio en la media (Cramer). Estas pruebas se pueden consultar 
en WMO (1971), y Machiwal y Jha (2012). Todas las pruebas citadas 
se aplicaron con un nivel de significancia (a) del 5 % y sus resultados 
se tienen en la Tabla 3; se indica con NH la pérdida de homogeneidad 
detectada. 


Tabla 3. Dos parámetros estadísticos y resultados de las siete 


pruebas de homogeneidad en las 22 estaciones hidrométricas de la 


Región Hidrológica No. 10 (Sinaloa), México. 


Pruebas de Ma 
o 
br sá Ba! el ha ll a hal 
al 


ps ei IA ad A ll ha 1 


ISO Ito hi sii IA li IO II 
a AAA AA 

Palo Dulce 0.305| -0.021 AA 
US ic lcd ad dio 
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Han ía 
H 


H|HAÍ]H 


al col E 1 


[ES H| HA NH] H 


Descripción de acrónimos: 


Cv = coeficiente de variación (S/X), adimensional. 
ri = coef. de correlación serial de orden uno, adimensional. 
CVN = cociente de Von Neumann. 


TA 


test de Anderson. 

PS = prueba de Spearman. 
TS = test de Sneyers. 

TB = test de Bartlett. 


PK = prueba de Kendall. 
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PC = prueba de Cramer. 


Se encontró que las siguientes seis estaciones hidrométricas 
muestran tendencia ascendente: Huites, Santa Cruz, Naranjo, 
Badiraguato, Zopilote y Pericos. Huites, Naranjo y Pericos además 
presentan persistencia, junto con las cinco siguientes: Guatenipa II, 
Jaina, La Huerta, Urique Il y Choix. Se observa que la persistencia está 
asociada con valores positivos del coeficiente de correlación serial de 


orden uno (ri) superiores a 0.20, excepto en Badiraguato. 


Prueba de discordancia 


Definidos los 22 registros por procesar (Tabla 2), se aplicaron las 
ecuaciones (2) a (9) para obtener los cocientes de momentos L, 
mostrados en la Tabla 4. 


Tabla 4. Cocientes de momentos L en las 22 estaciones 
hidrométricas de la Región Hidrológica No. 10 (Sinaloa), México. 


heal o MOM <A 
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AE EA 
2 San Francisco 0.23810 0.19358 0.12958 
3 0.22581 0.10956 0.17862 


NN Guatenipa II 0.23597 0.10004 0.11245 


Jaina 0.23033 0.17797 0.18571 


ESA Palo Dulce 0.17121 0.24359 0.04599 


Ixpalino 0.21747 0.02728 0.10958 


La Huerta 0.26480 0.27554 0.15224 


NN a an Dl 


AS A 
Vil Tamazula 0.18244 0.11958 0,14194 


14 | Guamúchil 0.36576 0.37388 0.25532 
15 | Choix 0.21305 0.18620 0.22018 


Badiraguato 0.32196 0.27868 0.20039 
El Quelite 0.48463 0.59616 0.56783 


18 | Zopilote 0.49516 0.33737 0.21933 
19 El Bledal 0.36091 0.22302 0.08793 
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21 0.45387 0.28892 0.16375 


En seguida, con base en los valores de la Tabla 4, y las 
ecuaciones (10) a (12), se desarrolló la prueba de discordancia, con la 
cual se obtuvieron los resultados expuestos en la Tabla 5. La prueba 
encuentra que los registros de Palo Dulce (Dé = 3.62) y El Quelite 
(Di7 = 4.12) exceden el valor crítico de la discordancia de 3.00 (Tabla 
1) y por lo tanto son discordantes. El registro de Pericos muestra un 
valor alto (D20 = 2.71) y por ello también fue eliminado. 


Tabla 5. Resultados de la prueba de discordancia (D;) en las 22 
estaciones hidrométricas de la Región Hidrológica No. 10 (Sinaloa), 
México. 

o AA Ñ 
A pales 
San Francisco 0.25 13 | Acatitán 0.27 
A A 
o 


Jaina 0.32 16 | Badiraguato 


Palo Dulce 3.62* 17 El Quelite 4,12* 


41 


Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua: ¡8 pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


o) y) Check for updates 
OPEN ACCESS q = ] 


Tecnología y 


CienciaszAgua 


A A 


Urique II 0.26 La Tina 1.06 


*excede el valor crítico. 


En relación con el cociente r,, se observa en la Tabla 4 que el 
valor de Palo Dulce es muy bajo, y los de El Quelite y Pericos muy altos, 
en comparación con sus valores circundantes por tamaño de cuenca, 
del orden de 0.22 y 0.33, respectivamente. Porotra parte, en la Tabla 
3 y con respecto al valor del Cv, se aprecia que la magnitud de Palo 
Dulce es muy baja, y las de El Quelite y Pericos bastante elevadas; en 


resumen, discordantes. 


Homogeneidad regional 1: índices de estacionalidad 


Para cada año de los 19 registros por procesar (Tabla 2) se detectó el 
mes con máximo escurrimiento. Tales datos se procesaron con base en 
las ecuaciones (13) a (16), y se obtuvieron los resultados que se han 


concentrado en la Tabla 6 y se tienen dibujados en la Figura 2. Se 
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observa una gran estacionalidad, pues en la mayoría de los registros 
su dirección media (9) ocurre en la segunda quincena de agosto y 


primera de septiembre, con ángulos que fluctúan entre los 230.9% en 
Bamícori (asterisco) y los 259% en Urique II (estrella negra). 
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Figura 2. Representación polar de los índices de estacionalidad en 
las 19 estaciones hidrométricas de la Región Hidrológica No. 10 


(Sinaloa), México. 


44 
JW 2022, Instituto Mexicano de Tecnología y ciencias del agua, ISSN 2007-2422, 
Tecnología del Agua. O pen Access bajo la licencia CC BY- 13(4), 01-74. DOI: 10.24850/j-tyca-2022-04-01 
NC-SA 4.0 (https ://creativecommons.org/licenses/by- 

nc-sa/4.0/) 


o) y) Check for updates 
OPEN ACCESS e . ps) 


Tecnología y 


CienciaszAgua 


Tabla 6. Índices de estacionalidad en las 19 estaciones hidrométricas 
de la Región Hidrológica No. 10 (Sinaloa), México. 


A] pen 1 [A 5996 A E dá] 3 cid 7819 


cio 4 Passa 8444 [dl ci ñ 7 Pisa 8476 
Francisco 


[an Cruz 242. [een aid 9055 al a pesa 3 dl 8857 


lua ld o Ia pe po _ 


WA pan 3 nad 5105 Ar dl 7 iaa 7626 


EME Hesse 3 hescal 8479 A scid 8 Meses 8572 


AM] 246.0 | 0.5544 El Bledal 232.5 | 0.8817 


Toahayana 255.0 0.6046 18 La Tina 245.1 | 0.7836 


[RO Urique II 259.0 0.5474 230.9 | 0.8121 
sii ic MEGA POCO AN ER) 


Respecto a la variabilidad, la menor dispersión se tiene en Santa 


Cruz (cuadro negro-blanco) y la mayor en Jaina (rombo negro). 
También se observa que seis estacionestienen Índice de estacionalidad 
del orden de 0.60 o menor, las cuales son Huites (cuadro negro); 
Guatenipa II (círculo negro-blanco); Jaina, La Huerta (rombo negro- 
blanco); Toahayana (pentágono negro), y Urique II. Estas seis 
estaciones tienen cuencas grandes mayores de 4 000 km?, 
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Con base en los resultados de la Tabla 6, mostrados en la Figura 
2, se concluye que los 19 registros procesados de VEA forman una 
región homogénea desde un punto de vista hidrológico. 


Homogeneidad regional 2: linealidad de los 
momentos ordinarios 


A cada uno de los 19 registros por procesar (Tabla 2), se les aplicó la 
Ecuación (17) para obtener sus momentos ordinarios de orden 1 a 4; 
a estos valores se les aplicó el logaritmo natural para obtener 
magnitudes más fáciles de dibujar. Tales resultados se muestran en la 
Tabla 7. 


Tabla 7. Logaritmos naturales de los momentos ordinarios (MOs) en 


las 19 estaciones hidrométricas de la Región Hidrológica No. 10 


(Sinaloa), México. 


E O 


8.370 16.883 25.917 34.250 
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San Francisco 7.812 15.788 23.908 32.149 


ASE 7.412 14.980 | 22.696 30.559 
HE Guatenipa II 7.355 14.867 22.506 30.251 


Jaina 7.224 14.610 22.138 29.787 


Ixpalino j 14.715 22.232 29.820 


Toahayana . 14.138 21.370 28.678 
Urique II ; 12.400 18.936 25,633 
10 Tamazula A 12.955 19.566 26.249 


11 Naranjo 5.308 11.264 17.656 24.332 


14 Choix 5.686 11.513 17.475 23.567 
Hill Bici ai id Bel 

16 Zopilote 

17 El Bledal 

18 La Tina 

19 Bamiícori 5.370 8.387 11.529 


Ordenada al origen (b) -3.723 -6.554 -8.837 -10.787 


Pendiente (m) 2.823 5.466 8.010 10.502 


Coef. correlación (r;,,) 0.964 0.965 0.964 0.963 
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En la Figura 3 se tiene una gráfica semilogarítmica, en cuyo eje 
de las ordenadas se han dibujado los momentos ordinarios de la Tabla 
7, y en las abscisas se tiene una escala logarítmica para las áreas de 
cuenca. Los valores de b y m de la Ecuación (18) se muestran en la 
columna respectiva de s en la Tabla 7. La gráfica interna de la Figura 
3 muestra la linealidad de las pendientes (m) y los valores de rx 
ratifican que las 19 estaciones hidrométricas analizadas forman una 
región homogénea hidrológicamente. 
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Pendiente (m) 


1 2 3 4 
Orden del momento ¿ o 


Logaritmos naturales de los momentos ordinarios (MOs) 


100 1000 10000 40000 
Área de cuenca (4), en km? 


Figura 3. Relación entre los momentos ordinarios y el área de 
cuenca en las 19 estaciones hidrométricas de la Región Hidrológica 
No. 10 (Sinaloa), México. 
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Homogeneidad regional 3: regresión lineal de los 
momentos f,. 


Ahora, a cada uno de los 19 registros por procesar (Tabla 2), se les 
aplicó la Ecuación (9) con r de 1 a 4, pues para r = 0 es la media y se 
tiene en la Tabla 2. Tales resultados se exponen en la Tabla 8. 


Tabla 8. Momentos de probabilidad ponderada (f,) en las 19 


estaciones hidrométricas de la Región Hidrológica No. 10 (Sinaloa), 


México. 


pa TE 640.703 | 1943.151] 1551.143] 1 id 210 
sra 1528.455 | 1135.915| 914.060 ro 370 


HG ca Noi sc ii ci 
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7 La Huerta 700.858 529.626 431.571 366.896 


EN 670.252 489.076 | 387.691 322.538 


A an 281.128 | 214.070 | 175.137 | 149.325 
366.742 | 265.611 | 210.191 | 174.859 


cl 149,155 120.729 102.637 90.002 


270.775 205.468 | 167.983 143.351 
87.910 69.388 58.374 50.916 


14 | Choix 178.797 131.614 105.574 88.851 
15 | Badiraguato 167.162 128.794 106.354 91.385 


16 | Zopilote 62.245 50.686 43.298 38.132 


El Bledal 31.054 24.060 19.812 16.934 


19 8.446 6.451 5.267 4.473 


Ordenada al origen (b) 7.4878 3.2501 2.0339 1.4745 
Pendiente (m) 0.6110 0.7352 0.7980 0.8362 


"y 


Coef. correlación (r;y) 0.9999 0.9999 0.9999 0.9999 


Los resultados de la aplicación de la Ecuación (19) se muestran 


en la columna respectiva de r de la Tabla 8, de manera que en la 


primera se tiene la relación lineal entre B1 y Bo, y en la última, la que 
existe entre Ba y Baz. 
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Con base en el valor del rxy (Ecuación (20)), se observa que tales 
relaciones son casi perfectas. Por lo anterior, este método también 
orienta a concluir que las 19 estaciones hidrométricas de la Tabla 8 


forman una región homogénea desde una perspectiva hidrológica. 


Eliminación de registros no homogéneos 


En la Tabla 3, los resultados de la aplicación de las pruebas estadísticas 
de homogeneidad local o individual mostraron que seis estaciones o 
registros procesados tienentendencia ascendente: Huites, Santa Cruz, 
Naranjo, Badiraguato, Zopilote y Pericos. Esta última ya fue eliminada 
en la prueba de discordancia. En los dos primeros también se detectó 
persistencia y en los cinco siguientes: Guatenipa II, Jaina, La Huerta, 
Urique II y Choix. Estas siete estaciones con persistencia se analizan 


en el inciso siguiente. 


Con base en las ecuaciones (21) a (27) se probó, en los seis 
registros con tendencia, si ésta era significativa; sólo en la estación 
Zopilote no lo fue. Debido a que la tendencia detectada es ascendente, 
se decidió eliminar tales registros, ya que sus datos tenderán a 
aumentar el escurrimiento anual y ello incrementará la incertidumbre 


asociada con su estimación mediante el análisis de frecuencias 
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regional. Por lo anterior y finalmente, ahora se procesan 15 registros 
de volumen escurrido anual, para definir su distribución de 


probabilidades regional. 


Respecto a la persistencia 


En los escurrimientos anuales, la persistencia es una propiedad 
estadística que se origina en la ocurrencia de años húmedos y secos; 
por ello, se compararon los periodos de registro procesado de las 
estaciones con valor mayor del coeficiente de correlación serial de 
orden uno (ri) de la Tabla 3: La Huerta (0.486), Urique II (0.381) y 
Guatenipa II (0.336), para buscar similitud ensus lapsos, que en efecto 
ocurre en los años de 1970 a 2000. Las otras dos estaciones que 
mostraron persistencia fueron Jaina (0.216) y Choix (0.297), con 
registros diferentes tanto en longitud como en ocurrencia. 


Khaliq, Ouarda, Ondo, Gachon y Bobée (2006) indican que la 
validez de las predicciones de crecientes obtenidas através del análisis 
probabilístico será cuestionable si las observaciones del registro por 
procesar son dependientes, es decir, tiene persistencia. Lo anterior es 
válido cuando se buscan crecientes o predicciones asociadas con muy 


bajas probabilidades de excedencia, pues se obtienen en el extremo 
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derecho de la distribución de probabilidades; caso que no ocurre al 
generar secuencias sintéticas de escurrimientos anuales, que abarcan 
toda la extensión de dominio del modelo probabilístico utilizado (Salas, 
1993; Cannarozzo et al., 2009). Debido al razonamiento anterior, no 


se consideró necesario eliminar las estaciones con persistencia. 


Selección de la FDP regional 


Todos los cálculos que se describen a continuación se desarrollan con 
datos adimensionales, pues los volúmenes escurridos anuales de cada 
uno de los 15 registros por procesar se han divido entre su valor medio 
o índice de escurrimientos (IE). A cada registro adimensional se le 
aplicaron las tres FDP recomendadas por Cannarozzo et a/. (2009), que 
son la Pearson tipo III (ecuaciones (28) a (37)): la Log-Normal de tres 
parámetros de ajuste (ecuaciones (38) a (45)), y la transformación 
potencial (ecuaciones (46) a (49)). Los resultados de la FDP que 
condujo al menor error estándar de ajuste (EEA, ecuación (50)) se han 
concentrado en la Tabla 9. Se observa que la FDP Pearson tipo III no 


alcanzó un mejor ajuste en ninguno de los 15 registros procesados. 
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Tabla 9. Mejores ajustes y valores estimados para los deciles 
indicados en las 15 estaciones hidrométricas de la Región Hidrológica 
No. 10 (Sinaloa), México. 


0.700 0.900 
San Francisco | TP2 | 0.060 | 0.528 | 0.733 | 0.919 | 1.151 | 1.591 
Jaina LN3 | 0.053| 0.532 | 0.743 | 0.927 | 1.150 | 1.558 
Ixpalino TP2 | 0.035| 0.516 | 0.789 | 0.988 | 1.194 | 1.502 


0.057| 0.514 | 0.695 | 0.873 | 1.117 | 1.663 


Toahayana TP2 | 0.059| 0.534 | 0.777 | 0.969 | 1.179 | 1.514 


Urique II OS 0.439 | 0.657 | 0.869 | 1.148 | 1.716 


Acatitán LN3 | 0.102| 0.412 | 0.662 | 0.890 | 1.175 | 1.722 
Zopilote LN3 | 0.166 | 0.115 | 0.402 | 0.713 | 1.164 | 2.199 
A A 
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En la Tabla 10 y Tabla 11 se muestra el detalle de resultados con 
las tres FDP en el registro con EEA más alto (El Bledal, peor ajuste) y 


con el más bajo (Tamazula, mejor ajuste). 


Tabla 10. Parámetros de ajuste y valores de los deciles en la 
estación El Bledal (la de mayor EEA) con las tres FDP citadas. 


PT3 LN3 TP2 
Parámetro 1 H = 1.000 €=0.857 p= 0.145 


Parámetro 2 o=0.677 a. = 0.585 Cs = 0.002 


EEA (adim.) 0.185 0.180 0.176 


P(X<x) V? /1E) V? J1E3 V?/1E) 


0.10 0.287 0.291 0.315 


0.20 0.438 0.449 0.445 


RETA 1.020 1.014 0.983 


0.70 1.217 1.204 1.181 


si dl > 
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Tabla 11. Parámetros de ajuste y valores de los deciles en la 
estación Tamazula (la de menor EEA) con las tres FDP citadas. 


PT3 LN3 TP2 
Parámetro 1 H = 1.000 ¿ = 0.961 p = 0.295 


Parámetro 2 o = 0.329 a =0.315 Cs = 0.001 


EEA (adim.) 0.035 0.035 0.034 


P(X<xX) VÍ JE v? /1E1 V? /1E) 


0.10 0.613 0.614 0.614 


0.20 0.719 0.721 0.722 


0.70 1,139 1.137 1.143 


Por último, en la Tabla 12 se exponen los resultados del método 


de las estaciones-años con las tres FDP utilizadas, es decir, 


conjuntando los 15 registros adimensionales, que suman 539 datos. 
Se concluye que la FDP Log-Normal conduce al menor EEA. 
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Tabla 12. Parámetros de ajuste y valores de los deciles en la Región 
Hidrológica No. 10 (Sinaloa), México, con las tres FDP citadas, 15 
estaciones hidrométricas y 539 datos. 


O 


. p 
AAA 
[AS 


Con base en los resultados de la Tabla 12, se define la curva de 


Pp LN3 T 
3 1) j 
l l 


T3 


P2 


crecimiento regional (Ecuación (1)), que permitirá generar secuencias 
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sintéticas de volumen escurrido anual adimensional en cualquier sitio 
o localidad de la Región Hidrológica No. 10 (Sinaloa), México; esta es: 


a(F)=0.896-22 [1 — exp(0.401-Z)] (52) 


0.401 


en la cual Z es la desviación normal estándarfunción de la probabilidad 
de no excedencia F. En la Figura 4 se han dibujado en el papel de 
probabilidad log-normal, los primeros 500 datos del volumen escurrido 
anual adimensional de 25 en 25 valores, después los siguientes 30 


datos de 5 en 5, y por último los cuatro valores finales, además de la 
FDP Log-Normal de tres parámetros de ajuste. 
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Periodos de retorno, en años 


10 20 50 100 1000 


MESA: 
OCRE 
ET 


il 2 5 10 20 30 40 50 60 70 95 98 99 99.5 9.9 
Probabilidad de no polo en pue de 


Figura 4. Ajuste de la distribución Log-Normal a los 539 volúmenes 
escurridos anuales adimensionales de la Región Hidrológica No. 10 
(Sinaloa), México. 
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Secuencias sintéticas de VEA adimensional 


Para la aplicación de la Ecuación (52), primero se debe definir el 
tamaño (ns) que tendrá la secuencia sintética generada de volumen 
escurrido anual adimensional. Después, haciendo uso de un generador 
eficiente de números aleatorios con distribución uniforme (Metcalfe, 
1997), se obtienen los ns números u; requeridos; tales valores se hacen 
igual a la probabilidad de no excedencia (F) y se aplica la Ecuación 
(30), para obtener los valores de Z necesarios, que se llevan a la 
Ecuación (52) para obtener las magnitudes q(F) buscadas. Con el 
proceso anterior se puede calcular cualquier número de secuencias 


sintéticas de V? J1EÍ (Ecuación (1)) igualmente probables de ocurrir. 


Escalamiento de las secuencias sintéticas 


Al dibujar en un papel logarítmico las dos columnas finales de la Tabla 
2, con el área de cuenca (4) en las abscisas y el volumen escurrido 
medio anual (VEMA) en las ordenadas, se define una línea recta que se 


muestra en la Figura 5 y cuya ecuación es: 
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IE! = VEMA= 0.015324- 41284914 (53) 


En la expresión anterior, el coeficiente de correlación resultó de 
0.966, con las 15 parejas de datos utilizadas. 


6000 


100 


Volumen Escurrido Medio Anual (VEMA), en Mm? 


4 Y " E 1 = 


100 1000 10000 


Área de cuenca (A), en km? 


Figura 5. Relación entre el volumen escurrido medio anual (VEMA) y 
el área de cuenca (4) en las 15 estaciones homogéneas de la Región 
Hidrológica No. 10 (Sinaloa), México. 
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La Ecuación (53) es válida en toda la Región Hidrológica No. 10 
(Sinaloa), y permite estimar el valor del índice de escurrimientos (1), 
con el cual se les asigna dimensión en Mm? a las secuencias sintéticas 
de VEA generadas con la Ecuación (52). 


Acopio de resultados principales 


De la Región Hidrológica No. 10 (Sinaloa), México, se procesaron 22 
registros de volumen escurrido anual (VEA). La prueba de discordancia 
detectó dos registros anómalos: El Quelite y Palo Dulce; además, se 


eliminó el de la estación Pericos, por mostrar discordancia elevada. 


En este estudio regional del VEA se evaluó la fecha de ocurrencia 
media (0) del volumen escurrido máximo mensual de cada año y su 
regularidad (7). En la Figura 2 se definen dos grupos de estaciones; las 
de regularidad baja con 0.60 o menos tienen cuencas grandes mayores 
de 4 000 km?. El resto de las estaciones tiene regularidad alta entre 


0.78 y 0.90. Sin embargo, por su fecha de ocurrencia media forman 
una región homogénea. 
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Las otras dos pruebas de homogeneidad regional, cuyos 
resultados se muestran en la Figura 3 y Tabla 8, ratifican que los 19 
registros procesados pertenecen a una región homogénea desde un 


punto de vista estadístico. 


Antes de iniciar la búsqueda de la FDP regional se eliminaron 
cuatro registros de VEA por mostrar tendencia ascendente 
significativa: Huites, Santa Cruz, Naranjo y Badiraguato. En la Tabla 2 
se observa que tales estaciones tienen áreas de cuenca diferentes, las 
dos primeras grandes y las otras pequeñas. Además, de la Figura 1 se 
deduce que por su ubicación no definen una subregión específica, en 
la cual se pudiera aceptar que tal componente determinística existe. 
Por lo anterior, se procedió a eliminarlas, ya que estarían influenciando, 


por exceso, las estimaciones de las predicciones. 


El procesamiento de los 15 registros de VEA restantes se llevó a 
cabo conjuntando tales registros para aplicar el método de las 
estaciones-años, pues con tales registros no puede definirse un periodo 
común de datos. Los datos anuales (VEA) se dividieron entre el 
volumen escurrido medio anual (VEMA) de cada registro, para ser 


procesados como valores adimensionales. 


Se probaron tres modelos probabilísticos: las distribuciones 
Pearson tipo III (PT3), la Log-Normal de tres parámetros (LN3) y la 
técnica estadística de la transformación potencial (TP2), que equivale 
a la FDP Normal. Se adoptó la distribución LN3 por conducir a un menor 
error estándar de ajuste. 
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Por último, se definió una ecuación logarítmica entre el VEMA 
como variable dependiente y el área de cuenca de la estación 
hidrométrica como variable independiente. Tal ecuación permite el 


“escalamiento” de las predicciones adimensionales obtenidas con la 
FDP regional. 


Comparación con otros estudios 


Los análisis de frecuencias regionales se han realizado de manera 
preponderante para las crecientes tanto en México como en el mundo. 
En cambio, el método de índice de escurrimientos (IE), descrito en 
este trabajo, no ha sido aplicado en nuestro país y a nivel mundial, la 
necesidad de estimar volúmenes escurridos anuales para el diseño 
hidrológico de los embalses ha sido resuelta mediante métodos que 
simulan el proceso lluvia-escurrimiento (Beven, 2001; Wagener, 
Wheater, 8 Gupta, 2004). 


La aplicación del análisis de frecuencias regional para estimar el 
volumen escurrido anual, en Sicilia, Italia, por Cannarozzo et al. 
(2009), sustentó el objetivo de este trabajo. Sin embargo, ambos 
estudios son completamente diferentes. En Sicilia se definieron tres 


regiones hidrológicas con base en la altitud y la lluvia media anual en 
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cada cuenca; su verificación se realizó a través de la prueba H de 
heterogeneidad (Hosking € Wallis, 1997; Campos-Aranda, 2010). 


Por otra parte, del análisis de frecuencias de crecientes realizado 
en Sicilia por Ferro y Porto (2006), se tomaron las pruebas de Gupta 
et al. (1994) y Kumar et al. (1994), para verificar la homogeneidad de 
la Región Hidrológica No. 10 (Sinaloa), México. 


Conclusiones 


Se ha desarrollado un análisis de frecuencias regional procesando 22 
registros de volumen escurrido anual (VEA) de la Región Hidrológica 
No. 10 (Sinaloa), México. Para seleccionarlos registros que forman una 
región homogénea, según un enfoque hidrológico, se aplicó la prueba 
de discordancia y se eliminaron tres registros. 


La verificación de la homogeneidad regional hidrológica de los 19 
registros restantes se realizó con base en tres pruebas estadísticas: 
(1) la de los índices de estacionalidad; (2) la linealidad de los 


momentos ordinarios, y (3) la regresión lineal de los momentos de 
probabilidad ponderada. 
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Con base en siete pruebas de homogeneidad local, se eliminaron 
otros cuatro registros que mostraron tendencia lineal ascendente. Por 
lo anterior, se procesaron finalmente 15 registros de VEA, con 
amplitudes de registro que variaron de 24 a 57 años, y sumaron 539 


datos. 


Para definir la curva de crecimiento regional [q(F)] se aplicaron 
tres modelos probabilísticos: (a) Pearson tipo III, (b) Log-Normal y (c) 
Normal mediante la transformación potencial. Los dos primeros 
ajustados con el método de los momentos L. Se encontró, a través del 
menor error estándar de ajuste, que la distribución Log-Normal 
representó mejor los 539 valores adimensionales de escurrimiento 
anual, y definió la ecuación siguiente, función de Z, que es la desviación 
normal estándar: q(F) = 0.896 — 0.497 - (1 — exp(0.401 - Z)/0.401 


Finalmente, para el escalamiento de las secuencias sintéticas 
generadas con la ecuación anterior, se estableció la relación empírica 
siguiente: VEMA = 0.015324 + 41284914. 


Para cuencas sin hidrometría, resulta imprescindible contar con 
un método que permita estimar registros de VEA de cualquier tamaño 
e igualmente probables de ocurrir que un registro histórico. Con la 
ventaja adicional para tales secuencias sintéticas de tener por su 
naturaleza aleatoria ausencia de persistencia y con ello conducir a un 


dimensionamiento de los embalses más confiable. 


Por lo anterior, se recomienda la aplicación sistemática del 


método del índice de escurrimiento, según el procedimiento expuesto, 
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en otras regiones hidrológicas del país para encontrar las ecuaciones 
regionales correspondientes de q(F) y VEMA = f(:). 
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